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I gauge-invariant thermal renormalisation group. 
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An important goal in thermal field theory is the computation of physical quantities like 
the pressure or screening masses of quantum fields. This is of interest for on-going and 
future experiments at RHIC, LHC and the FAIR facility, which are sensitive to the quark- 
gluon-plasma in the core region of heavy ion collisions. Recent RHIC data even suggests 
that a thermal quark gluon plasma is strongly coupled, see [1] . On the theoretical side, one 
observes that standard thermal perturbation theory is plagued by infrared divergences due to 
massless excitations, which at least requires resummations. Resummed thermal perturbation 
theory often displays a weak convergence behaviour even deeply in its domain of validity, 
see, e.g. [2-4]. In gauge theories, the situation is additionally complicated by the magnetic 
sector. In our opinion, this asks for a method able to incorporate both the weakly and the 
strongly coupled regimes. Moreover, such a method would benefit from a clear separation 
of quantum fluctuations and thermal effects. 

In the past decade much progress in the description of quantum fluctuations has been 
achieved within the functional renormalisation group, for reviews see [5-11]. The functional 
RG, applicable both at weak and strong coupling, has also been used for the computation of 
thermal effects, e.g. [6-8, 12-26]. A direct implementation of thermal fluctuations has been 
put forward within the real time formalism [14-16] and the imaginary time formalism [7]. 
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In the present paper, we extend and detail the proposal [7]. We construct thermal 
flows that are well-defined both in the infrared and in the ultraviolet, and provide ideal 
starting points for numerical and analytical studies. In scalar theories, we work out the 
scheme-independent part of the integrated flow, which contains the one loop thermally 
resummed perturbation theory. We equally provide thermal flows within a thermal derivative 
expansion, and reproduce explicitly the leading order results for the pressure and the thermal 
self energy. 

For gauge theories, and for specific momentum cutoff and gauges, the thermal fiow even 
respects gauge invariance for arbitrary scale k. This follows from two observations. A 
mass-like momentum cutoff - even though it is not applicable for integrating-out quantum 
fiuctuations [7] - is a viable cutoff for thermal fluctuations. The second observation concerns 
Wilsonian flows for gauge theories in axial gauges [27, 28]. Gauge invariance for physical 
Green functions is controlled via modified Ward-Takahashi identities. In an axial gauge they 
reduce to the standard Ward-Takahashi identities for any scale A;, if a mass-like regulator is 
employed. Although of limited use in the generic case, this is precisely the missing piece to 
construct a gauge invariant flow for thermal fluctuations. 

The outhne of the paper is as follows. In Sect. II, we review functional flows in the 
imaginary time formalism. In Sect. Ill, we discuss our projection method for thermal fluc- 
tuations. Scheme independence of the full one-loop thermal resummation is worked out in 
Sect. IV. Explicit thermal flows within a thermal derivative expansion and results for the 
thermal pressure and the self-energy are obtained in Sect. V and Sect. VI. The extension 
to thermal gauge theories is detailed in Sect. VII, including a gauge-invariant thermal flow 
and results for the thermal mass and pressure. Our conclusions are contained in Sect. VIII. 



To start with, we discuss the standard Wilsonian approach to quantum fleld theories. For 
simplicity, we focus on a bosonic degree of freedom. A generalisation to fermions and gauge 
flelds is straightforward. All the physically relevant information can be obtained from the 
(regularised) partition function. It reads 



In {d + 1) dimensions and using the imaginary time formalism at temperature T, the trace 
stands for 



II. FUNCTIONAL FLOW EQUATIONS 




(1) 




(2) 



and the implicit replacements go ~^ 2TinT for bosonic flelds are understood, n labelling the 
Matsubara frequencies and J stands for the corresponding sources. The term Sk,T[4>\ = 
S + ASk,T contains the (gauge-flxed) classical action S[(p] and a quadratic regulator term 
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A5'fe,r[0], given by 



ASkM = H-q) Mq') Hq) ■ (3) 



Eq. (3) introduces a coarse-graining via the operator Rk{q)- The flow of (1) related to an 
infinitesimal change oi t — Ink/ A (with A being some fixed UV scale) is 

^ ^w,,r[J] = 1 y p0 Tr (^0 ^ 0^ e-^^--['^l+^ ■'f' . (4) 

Performing a Legendre transformation leads to the coarse-grained effective action rfc^r[0], 

rk.T[<P] = Tr J0 - Wk.riJ] - A^fc,r[0] - C^.t , (5) 

Note that the constant Ck is usually not mentioned when one is interested in field indepen- 
dent quantities. However as we shall see later in the discussion of the thermal pressure we 
have to take it into account. It is straightforward to obtain the flow equation for rfcr[0] by 
using (5): 



dtrk,T[<l>] = ^Tr 



- dtCk,T , (6a) 



with 



denoting the full (fleld-dependent) regularised propagator of 4>. Let us be more speciflc 
about the regulator function Rk{q)- The general requirements are 

(i) it has a non-vanishing limit for 0, typically Rk — ^ for bosons, and Rk ^ k for 

fermions. This precisely ensures the IR flniteness of the propagator at non-vanishing 
k even for vanishing momentum q. 

(ii) it vanishes in the limit k ^ 0, and for q"^ ^ k"^. The latter condition ensures that large 
momentum fluctuation have efficiently been integrated-out whereas the first condition 
guarantees that any dependence on Rk drops out in the limit A; — > 0. 

(iii) Rk diverges like for bosons, and like — >■ A for fermions, when k ^ oo (or 
k ^ A with A being some UV scale much larger than the relevant physical scales). 
Thus, the saddle point approximation to (1) becomes exact and r^^A reduces to the 
(gauge-fixed) classical action S. 

These conditions guarantee that has the hmits 

hm TkM = S[^] (7) 

re— >oo 

iimrfe,T[0] = rT[0] . (8) 
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For any given scale k the main contributions to the running of T^ t in (6) come from momenta 
about ^ A:^, if the regulator function obeys the conditions This is so because 

dtRk is peaked about ^ k^, and sufficiently suppressed elsewhere. The physics behind 
this is that a change of ^ due to a further coarse graining (i.e. the integrating-out of 
a thin momentum shell about k) is dominated by the fluctuations with momenta about 
k. Contributions from fluctuations with momenta much smaller/larger than k should be 
neghgible. The flow equation (6) connects the classical action S[(j)] with the full quantum 
effective action Fr[0] at temperature T. The full quantum effective action is hence defined 
by the fiow equation and an initial action at some UV scale. Note that the limits (7) and 
(8) are strictly speaking only valid with a suitable choice of C^. 



III. QUANTUM VS. THERMAL FLUCTUATIONS 



At finite temperature the fiow (6) constitutes the integrating-out of quantum as well as 
thermal fluctuations. Eq. (6) can be projected on thermal fluctuations: instead of computing 
the flow for Ffe^^f^] as in (6), we propose to study the flow for the difference 

f fc,T[0] = - rfe,o[0] (9) 

between the effective action at vanishing and non- vanishing temperature [7]. Eq. (9) is 
evaluated for flelds periodic in time.^ It entails the thermal effects in the effective action as 
we have effectively removed the quantum fluctuations. The flow of (9) reads 

dthA^P] = ^ / tSi E ^i^M dtRk - / ^ G,,o[0] dtRk] - dtCk,T , (10) 

where Ck,T = Ck,T — Ck,o- The fiow (10) vanishes at T = 0, and hence is a thermal fiow. 
It can be used within a two-step procedure: first, one computes the fiow (6) at vanishing 
temperature within a truncation that is adapted to the zero temperature theory. Then the 
result is used as an input for the thermal flow (10). We emphasise that the truncation at 
flnite temperature may differ even qualitatively from that at T = 0. 

The thermal flow (10) has another important advantage that is particularly interesting 
for its application to gauge theories. We can relax the conditions (i) — (iii) on the regulator. 
The aim is to flnd a reliable, but still sufficiently simple and manageable formulation of 
the fiow equation. The first observation is that the flow equations (6) and (10) indeed are 
simplified for a mass-like regulator given by 

bosons: Rk{q) ~ k^ 
fermions: Rk{q) ~ 



For Thfi, this means that Greens functions F^ q are multiphed with periodic fields, and integrated over 
the compact time interval. 
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or variants of it.^ The key characteristic of Rk in (11) is that it does not depend on momenta. 

Formally speaking, (11) is a viable IR regulator in the sense of condition (i), and it 
allows as well to reach the UV initial condition, due to condition {iii). However, the choice 
(11) violates condition {ii), which is one of the basic requirements for a Wilsonian cutoff. 
Indeed, the operator dtRk appearing in (6) is neither peaked about ^ A;^, nor does it 
lead to a sufficient suppression of high momentum modes. The flow equation (6) would then 
receive contributions from the high momentum region for any value of k. An immediate 
consequence of this is that an additional UV regularisation is required, as the flow equation 
(6) is no longer well-deflned for large loop momenta. Stated differently, on might say that 
a mass term regulator leads to a break-down of the Wilsonian picture, since it is no longer 
related to an integrating-out of momentum degrees of freedom. Rather, it corresponds to a 
flow within the space of massive theories. 

Apart from these more formal objections one should mention that numerical solutions of 
(6) are much more involved than for regulators satisfying {ii). At every iterative step a d- 
dimensional momentum integral has to be performed in (6) over a non-trivial function which 
is not strongly peaked about some momentum region. This is, numerically, a quite tedious 
problem. For this reason most of the sophisticated numerical investigations are based on 
non-local regulator functions (like the sharp cutoff, exponential or algebraic ones). It has 
also been observed that approximate solutions to the flow equation (expansions in powers 
of the field, derivative expansions) show a rather poor convergence behaviour, when (11) is 
used. In gauge theories, this has been seen from explicit perturbative computations [29-31]. 
Consequently at T = it is not advisable to directly resort to a mass-like regulator. 

While (11) is not viable for quantum fluctuations being present in (6), we will now 
argue that it is perfectly viable for thermal ones in (10). The main point is that the large 
momentum fluctuations (not sufficiently controlled by (11), introducing UV divergences to 
the flow equation) have nothing to do with the heat bath. Therefore, subtracting the zero 
temperature quantities will render the flow equation (10) flnite and well-deflned, even in the 
presence of a mass-like regulator Rf. — k'^. Then the flow (10) boils down to 



Here, the momentum-independent regulator dtRk ~ Rk acts only as a multiplicative constant 
because of (11). In this case, the suppression of large momenta does not originate from dfRk, 
but from the cancellation between the propagator terms. For large internal momenta, the 
Matsubara sum can be replaced by an integral, thereby cancelling the T — contribution. 
Therefore, one may read (12) as a Wilsonian flow for thermal fluctuations: At the starting 
point = A (A being some large UV scale) all fluctuations are suppressed and (12) vanishes. 
For any k < A, the flow of rfc^r[0] would receive contributions for all momenta. In contrast, 
the difference f ^.^[0] is sensitive only to thermal fluctuations, which are peaked in the 



^ Sometimes it is convenient to multiply the mass term with a (momentum-independent) wave function 
renormalisation, Rk = Zk fc^ (and analogous for fermions). The same reasoning applies. 




27r 



(12) 
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infrared region and naturally decay in the UV region. It follows that the integrand in (12) 
is peaked about ^ k"^. In other words, condition (ii) is effectively guaranteed even in 
the case of a mass-like regulator by the very nature of the temperature fluctuations. This 
amounts to the fact that the mass-like regulator seems to be a reasonable choice even for 
numerical applications in thermal field theories. 

It is worth mentioning that the cancellation of UV divergences in (12) is reminiscent of 
the BPHZ-procedure. There, the subtraction of possibly divergent terms takes place on the 
level of the integrand rather than on the level of the regularised full expressions. 

Let us comment on the initial condition to (12) [and (10)]. In contrast to the flow (6) 
with the hmits (7) and (8), the flow equation (12) [resp. (10)] has the limits 

lim = (13a) 

iimffc,T[0] = rT[0] - ro[(/)] . (13b) 

The boundary condition (13a) looks rather simple. The flow equation (12) needs in addition 
the knowledge of the massive T = quantum theories. This point is qualitatively shared 
by the proposals [14, 38]. It seems likely to find a good approximation for the issues under 
investigation, since (12) is eventually projecting-out thermal fiuctuations. Those should 
not be too sensitive to the details of the quantum effective action at T = 0. Moreover we 
deal with a situation were the original fields are still sensible degrees of freedom. Thus, a 
perturbatively resummed quantum effective action should be a good starting point. Here we 
are actually taking advantage of the fact that for a mass-like regulator the flow is describing 
a path in the set of massive vector boson theories rather than a Wilsonian integrating-out. 



IV. THERMAL RESUMMATIONS 



The present formulation naturally incorporates thermal resummations within general 

truncation schemes. Moreover, an important, direct, consequence of the flow (10) is scheme 
independence of lowest order resummed perturbation theory: at this order the effective 
action at /c = does neither depend on the regulator nor on the truncation. This is proven 
below, and put to work within the derivative expansion in the next section. We discuss a 
scalar 0^-theory with a massless (at T = 0), neutral scalar field and coupling A, 



(14) 



where = j^^^ dr J d^x. We are interested in the A-dependence of thermal corrections. 
It is well-known that naive perturbation theory breaks down beyond one-loop, and we are 
left with an expansion in A^/^ rather than in A. For the proof of scheme independence we 
rewrite the fiow (10) as 



dtTk,! 



iTr 
2 



r{2) 

k,T 



Rk 



<2) 
0,T 



+ Rk 



9A + i^p5y^ 

0,T 



1 1 

dfRk — — Tr — T-r — 

+ Ri. 2 r^^^ + 

-t rLk i j^ Q -h 



Rk 



-dtRk 



(15) 
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where we have dropped the constant contribution. The traces Tr in (15) have to be taken 
at temperature T (first terms), and T = (last term). For computing the one loop effective 
action from (15) we have to insert the classical action on the rhs of (15). Then the first 
term on the rhs of (15) vanishes. The remaining terms are total t-derivatives and can be 
integrated trivially. After a reordering we end up with 



pl— loop 
^ k,T 



"-"cl.T 



o„i rp -|- n,ft 



loop 



-Tr In 



/pl— loop pl~l 
~ V- k,0 ~ ^ A,0 



loops 



(16) 



The second bracket in (16) vanishes identically due to (6). The A, T-dependent terms in (16) 
just arrange for the appropriate renormalisation; indeed their sum has to be A-independent. 
This results in 



pl— loop 
^ k,T 



\ ^ ct>{x) {d' + n;7^-p) 0(x) + A ^ <t>\x) + o{\') , 



(17) 



where 11^ 



nin 



a'^{l + 0{k/T)) is the momentum-independent self-energy at one loop, and 
oc A is the standard thermal mass at k = 0. Its scheme independence follows from (16). 



At k 



pl— loop 
^ 0,T 



we are left with the renormalised one-loop effective action at finite temperature 
= rQ^'°°^-|-ArQ Now we re-insert r^^"°''', (17) on the rhs of (15), and project the 
flow (15) on the leading resummed term of order A^/^. This term only receives contributions 
from the vanishing Matsubara frequency n — 0. Due to the subtraction the first term on 
the rhs of (15) has an additional (dimensionless) power of Xk/T. If integrating the zero 
frequency flow from A — > cxo to A; = 0, we can rescale all expressions with the one loop 
self energy H^t = ""^t- The first term in (15) is of order A^ due to the additional power 
Xk/rriT oc A^^^, whereas the second term receives A^'^^T^-contributions from the thermal 
trace at /c = 0. This term is a total i-derivative and we arrive at 

r(2) 



1. 



-Tr In 



0,T 



^(2) 



+ rA,T + e'(A'), 



(18) 



where the 0(A^)-terms also include additional subtractions. Eq. (18) includes the full 1- 
loop resummed perturbation theory, but also applies to more general truncations of ro,T- 
We emphasise that (18) is valid for general regulators and truncations and hence is scheme 
independent. Due to the scheme independence of the one loop effective action Tq^t, (18) 
scheme-independently contains the full bubble resummation that leads to the lowest order 
in the thermal resummation. We conclude that the flow (10) leads to the full thermal re- 
summation at lowest order for any regulator and any truncation. This includes in particular 
the thermal mass contribution at order A^^^ as well as the pressure contributions at order A 
and A^/^. In scalar theories this is a direct consequence of the fact that these terms are only 
related to bubble diagrams. 



V. THERMAL DERIVATIVE EXPANSION 



It is useful to discuss the above results in a very simple truncation, the leading order in 
the derivative expansion at finite temperature. In this approximation, the effective action 
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for a real scalar field reads 



(19) 



with the effective potential Vk,T{<P)- We follow the reasoning of Sect. Ill and consider the 
flow for the thermal difference Uk,T{4>) = Vk,T{4>) ~ ^fc,o(0), explicitly given by 



dtUka 



(2^ 



dtRk 



dqo 



dtRk 



(20) 



We have dropped a field-independent constant, and q^ — q^ + q^- For lowest order resummed 
perturbation theory we can concentrate on the vanishing Matsubara frequency. The flow 
(20) then reads 



dtUk,i 



-T 



cPq 



dfR 



(27r)3q2 + i?, + T/^';^ 



+ C>(A2) = -T 



cPq 



dfR 



(27r)3q2 + i?, + T/o';^ 



+ 0{X'), (21) 



and (18) follows. 



At finite temperature a regularisation of the theory is already achieved for cutoffs depend- 
ing only on spatial momenta, Rk = Rk{<i)- These cutoffs are well-adapted to the present 
situation: then the thermal flows (10) successively sum up thermal fluctuation at momenta 

^ k"^, and the thermal properties of the theory are unchanged. Technically this has the 
benefit that the Matsubara sum is unchanged. In the present case this allows us to perform 
the Matsubara sum and the go-integration in (20) analytically. 



1 



^k,T l^kfi 



+ 



<^k,T 



(22) 



where u;k,T = (q^ + -Rfe(q) + V^'j,) 
bution 



1/2 



We note the appearance of the Bose-Einstein distri- 



n{u!) 



(23) 



in the flow, a direct consequence of the thermal sum. Every single term on the rhs in (22) 
remains well-deflned even for large spatial momenta, due to dtR. Along the flow, (22) uses 
the zero temperature running for \4,o- While the last term in (22) is clearly of a thermal 
origin, the flrst two terms display a non-trivial cancellation of zero temperature and thermal 
contributions. 

In the high temperature limit, the Bose-Einstein distribution is strongly enhanced, with 
n{u!) — > T/cu. Consequently, the flow (22) is parametrically dominated by the last term 
for T — > oo, and flxed k. In this limit, we introduce rescaled flelds as (p — (f)/VT and the 
rescaled potential as T4(</?) = Uk{(t>)/T. Then 



dtVk^ 



1 



cP'q 



dtR{q) 



2 J (27r)3q2 + i?,(q) + l-," 



(24) 
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All temperature dependence has disappeared and we are left with a standard flow in three 
dimensions. We emphasise that (24) equals the zero mode flow in (21), as all higher modes 
decouple due to their masses 27rT oo. We also point out that (24), unlike (22), is closed 
since contributions from u;^ o are absent in this limit. 

In the low temperature limit T ^ 0, the thermal flow (22) vanishes identically. At flnite 
temperature, and in a one-loop approximation, Vk^r is substituted by the classical potential 
V^i and the first two terms in (22) vanish. Then the thermal corrections reduce to a total 
^-derivative, and we obtain the cutoff-independent result 



(25) 



upon integration. This equals (16) in the present truncation. From (25), we derive the 
free pressure for a single bosonic degree of freedom P[T] = 7r^T*^/90 (where P = —V). At 



(2) 

this order, the self-energy correction to the propagator r^^(q 
independent and given by a thermal mass no,T(q^ = 0) 



n^y°°'^(q) is momentum- 



rrin 



AT^/24. We have proven in 
Sect. IV that the next-to- leading order of the self-energy is also scheme independent. This is 
most directly seen by studying (21) for the self-energy. The n — Matsubara mode generates 
a A^/^-term from the leading order thermal mass. The scheme independent coefficient follows 
from integrating the total t-derivative in (21), leading to Smj. — —AT 771^/(871). Hence, 



no,r — 



24 



1-3 



A 



247r2 



1/2 



+ 0{X) 



(26) 



The thermal pressure can be obtained by re-inserting the self-energy Ilk,T into the flow (20): 
r^ .p[0 = 0] = (27rnT)^ -|- -|- Ilk,T- This produces /c-dependent bubble diagrams with 
regulator insertions. The lowest order (two loop vacuum bubbles) are total ^-derivatives, 
and the field-independent part ro,T[0] reads 



90 



8 V 247r2 



15 

T 



247r2 



3/2 



+ O(A') 



(27) 



In the above discussion, we have focused on scheme independent aspects of thermal flows. 
Now we specify explicit momentum cutoffs, which should be useful for numerical or analytical 
studies beyond the present order. Following the reasoning of Sect. Ill, we consider first a 
mass term momentum cutoff Rk{(i) = k"^. Then (22) turns into 



Wk,T 



(2^ 



1 



.^k,T 



coth 



1 



(28) 



where u!k,T = (q^ + k'^ + Vk^TV^"^ depends on spatial momenta q. At fixed momentum cutoff 
k, we note that both terms in the integrand diverge quadratically in the UV, while the sum 
remains well-defined even for large loop momenta. Furthermore, terms sensitive to V^'j^, V^'q 
or their difference are suppressed by additional powers in 1/q^ and are, hence, finite. This 
is an explicit example of the cancellation emphasised above. 



10 



VI. OPTIMISED THERMAL FLOWS 

The physical content of a given truncation and its stabihty and convergence towards the 
physical theory is enhanced by suitably optimised choices for the Wilsonian cutoff [11, 32, 33]. 
Here, the thermal flows (22) and (28) still require a spatial momentum integration due to 
the non-trivial q-dependence of the integrand. We take advantage of a simple optimised 
momentum cutoff [34] 

i?,(q) = {e - <i')9{e - q^) , (29) 

which is known to display remarkable stability properties. This regulator cuts off the propa- 
gating momentum modes in (22) homogeneously, since q^+i?(q^) = k'^ becomes momentum- 
independent for all infrared modes with < k^. Furthermore, the remaining loop integra- 
tion in (22) is performed analytically, which is a crucial advantage in view of numerical 
studies. The result is 



k 

dtUk,T 



\ I k k \ k , . 

H n{uJk,T) 



2 \<^ik,T <^ik,0/ '^k,T 



(30) 



1 /9 

where cuk^r = {k^ + ^^'t) ■ The truncated flow (30) clearly displays the thermal flow 
structure: the last term is proportional to the thermal distribution n{Lo) which vanishes for 

T — > 0, and requires no UV rcnormalisation. The first term has the subtraction at T = 
which removes the quantum fluctuation. It is reminiscent of the UV-subtraction in thermal 
perturbation theory. Here, however, the subtraction is not necessary to get a finite flow. 
Without the subtraction, the flow (30) has been provided as a proper-time flow in [23]. As 
shown in [35], such a proper-time flow derives from (29) to leading order in the derivative 
expansion. Beyond this order, it carries inherent approximations discussed in [35]. 

Within the flow (30) the study of the scheme independent leading order is most convenient 
due to its analytic structure. We solve (30) iteratively, starting at Ua^t — Va,t — V^,o = 
and vanishing mass at T = 0. To leading order in A, corrections in Uk,T and its derivatives 
originate solely from the thermal (last) term in (30). Performing also the /c-integration we 
find the free pressure P — 7r^r^/90. For the running thermal mass, we find 



127r2 



n'{uJk,T) n{uJk,T) \ . 1 / Afc,T A 




\fe,0 

, .2 , ,3 / 9 I , ,3 ,T3~ 



(31) 



To leading order in the coupling, we have \k;r — ^kfl — ^-nd only the first term in (31) 
contributes. To this order, the running thermal mass reads 

A 

^M-ioop = dx [xn{x) - x^ n'{x)\ + m\^_^^^^ . (32) 

The integral can be performed in terms of poly-logarithms. For A — > cxd we have m\ ^ — > 0, 
and at A; = we find m\ — rnQ^^_-^^^^ — AT'^/24. While the /c-running in (32) is scheme- 
dependent, the endpoint is not, and coincides with the standard expression for the one- loop 
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thermal correction. Beyond the leading order, the running thermal masses have to be taken 
into account. In lowest order resummed perturbation theory the Matsubara zero mode leads 
to a A^/^-correction to the self-energy. This term is produced in the flow from scales 

A; < 27rr , (33) 

also assuming A ^ 1. In this region the flow effectively reduces to the high temperature 
hmit (24) with n{uj) — ^ T/uj, as T is larger than all other scales. We derive from (30) 



9tUk,T = ^ ( ^ + 0{k/u;,,T) 1 , (34) 



leading to 



d,mi^ = -X,,T—, ( ^ + O(fcMV) I • (35) 



for the mass. At scales (33) the running of the one loop mass is sub-leading and we can 
substitute rnl -^_i^^^ — > in (34), (35). Moreover we are only interested in the A^/^-term 
777.3^2- For dimensional reasons this term is proportional to Xk^TTrriT — XTuit + O(A^). 
After rescaling k with rriT we arrive at the integrated flow 



^3/2 



A^r^477A.-l)--A^> (36) 



which leads to (26). The subtraction of 1 removes the infrared one loop contribution. We 
proceed with the A^/^-contribution P3/2 to the thermal pressure P, which follows directly 
from (34), 



Tml r „ f k^ ,0 A Tml 



The subtraction oi k^ — 1 removes the tree level and one loop infrared contributions. Eq. (37) 
provides the A^/^-order in (27). 

In Sect. IV we have shown that the results (36), (37) are scheme- independent. Therefore 
it should be possible to turn the corresponding flows into total i-derivatives. We exemplify 
this structure with an independent derivation of the thermal pressure up to order A^/^. The 
flow dtPk,T of the pressure is given by (minus) (30) at vanishing field. It depends on the 
running mass via ujk^^{(f)) — {k"^ + m\ j, -\- |A0^)^/^ at = 0. Note that we keep a classical 
A for technical reasons. The pressure does not depend on this choice. Now we rescale the 
mass in uj^^^t with a parameter a, m\ rp —>■ am\ rp. With this modification we get 

P{a) = Pfree + aPi + a=^/'P3/2 + 0{\^) , (38) 

since mg^g does not contribute to P3/2, and the A-dependence only originates in mfe^i_ioop- 
Now we take the a-derivative, 

da dtPk,T = -mlT\dl dtUk,TU=o + C(A') ^ -j 9t^fe,T + , (39) 
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m'iTdtmlT = T^dtiml^T? ■ (40) 



where we have used that d^oJu^T = '^'i^T^lj, ^k,T at vanishing field. Note also that dtmi rp 
depends on A^^- — A only at the order (9(A^), as already indicated below (35). Taking 
derivatives as in (39) is a standard computational trick used in thermal perturbation theory 
where mass-derivatives help to analytically perform thermal sums and momentum integrals. 
Eq. (39) carries an a-dependence since dtml j.{a) is a function of a via u!k,T- At a = 1 the 
flow (39) can be expressed in terms of a total ^-derivative, 

1 
2 

Eq. (40) relates to the t-derivative of the two loop vacuum bubble, and contains the full 
contribution to the pressure of order A. This reduction leads to the diagrammatic structure 
of resummed perturbation theory: the flow can be written as a t-derivative of the infrared 
regularised resummed diagrams. Now we restitute the a-dependence in (40) with the help 
of (38). By taking the a-derivative of (38) we conclude that the one loop term in (40) 
qis a-independent, and the resummed contributions of order A^/^ are proportional to a^^^. 
Performing the ^-integration we arrive at 

d^Po,T = -^m^ + ^cV^ml + O(A^) = (^1 - ^c^^^) + 0(A^) , (41) 

where we have used (26) for substituting the A-dependence with the appropriate powers of 
rriT- We integrate (41) over a from to 1, and arrive at 



0,T 



^T--^ +0{y), (42) 



90 48 V TT T 
where we have used that -Po,rU=o is the tree-level pressure. Eq. (42) agrees with (27). 



VII. THERMAL GAUGE THEORIES 



We now turn to thermal flows for gauge theories. We aim at deriving well-deflned flows 
for gauge theories in terms of the fundamental degrees of freedom, the gluons. The inclusion 
of matter fields is discussed elsewhere. For specific regulators and gauges this flow is gauge 
invariant. The classical Yang-Mills action is given by 



SA[A] = ^j^d^xtTF\A), 



(43) 



with field strength F^j^^ = d^A^ — d^Ay -\- A„], and tntH^ — —5°-^/2 in the fundamental 
representation. The regulator term for gauge fields is given by 



^Su[A] = -TY, 



(44) 



where the imphcit replacement go 27mT is understood. With the regularisation (44) and 
using the definition (9), the fiow (10) reads 



(2^ 



UfJ, 

ba 



dtCk,T (45) 
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In the presence of matter fields and corresponding regulators additional terms pro- 
portional to dfRcj) will be present in (45). For mass-like regulators the flow (45) reduces 
to 



dtTk,T[A] = fc2 y tr TJ2 Gk,T[A] dtRk -J^ G,,t[A] 



k,T 



(46) 



where the trace tr sums over Lorentz and gauge group indices. So far we have not specified 
the gauge fixing procedure. In general, gauge fixing leads to additional ghost fields C as 
well as a rcgularisations Re- Since the thermal bath singles out a preferred rest frame, we 
have an additional Lorentz vector at our disposal. Therefore, it is natural to employ 
an axial gauge fixing. An axial gauge has the further advantage that ghost fields decouple 
completely from the theory, and possible Gribov copies are absent. We employ 



cfiq 



1 



4 



M f 2 "'i^'^v ■ 



(47) 



In [7, 27, 28] we discussed the various aspects of an axial gauge fixing (47). In particular 
we showed that the spurious propagator singularities of perturbation theory are naturally 
absent in a Wilsonian approach. Furthermore, the gauge fixing parameter ^ with mass 
dimension —2 has a non-perturbative fixed point at = 0. This singles out the nA = 
gauges and tremendously simplifies the problem of gauge invariance, because it allows for a 
momentum independent choice of ^. 

Gauge invariance for physical Green functions corresponds to the requirement of a mod- 
ified Ward Identity (mWI) to hold. Gauge transformations on the fields are generated by 



5a with action 



SaA 



D{A)a. 



(48) 



For momentum independent gauge fixing parameter ^, the mWI in the presence of the cutoff 
term (44) reads 



(49) 



The two terms on the r.h.s. are remnants from the gauge fixing and the coarse-graining, 
respectively. The mWI (49) is compatible with (6) [11, 27, 28, 36, 37], i.e. a solution to 
(49) at some scale k — h. remains a solution for A; < A if F^^T-f^] is integrated according to 
the flow equation. In particular, the terms proportional to R^ vanish for k 0, thereby 
ensuring gauge invariance for physical Green functions. The mWI related to the thermal 
flow (45) follows from (49) as 



SaTkAA] 



(2^ 



tr 



2n 



[a, Rk] GkAA] -TJ2[<^^ Rk] Gk,T[A] 



(50) 



The compatibihty of (50) with (45) is a direct consequence of the compatibihty of (49) with 
(6). The linear term related to the gauge flxing [the flrst term in (49)] has cancelled out, since. 
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as emphasized in Sect. Ill, we are only looking at fields (f) at temperature T and corresponding 
gauge transformations. This also implies - up to modifications for topologically non-trivial 
configurations - that a has to be periodic. Apart from this simplification, the same reasoning 
as for (49) above applies. 

With a mass-like regulator, however, we can go a step further. For a regulator as in (11), 
the right hand side in (50) vanishes since then [a, Rk] = 0. All coarse-graining dependence 
of (50) drops out for arbitrary scale k, and not only in the limit k ^ 0. This is an immediate 
consequence of Rk being momentum independent and the axial gauge fixing [27-29], and 
reduces (49) to the standard Ward Identity in the presence of an axial gauge fixing. It 
follows, that 

6aTk,T[A] = (51) 

and we end up with the statement that (12) corresponds to a gauge invariant thermal 
Wilsonian renormalisation group for Fk^Ti^'], valid at any scale k. 

The above flow can be used for computing the self-energy and the thermal pressure for 
one-loop resummed perturbation theory. Following the computations in Sect. IV and V we 
arrive at the well-known result for the thermal mass 



rr^l.on = l9'T'N. (52) 



The thermal mass (52) can be re-inserted into the flow. However, full resummations require 
also the one-loop running of propagator and classical vertices. This is straightforward but 
tedious, and the results shall be presented elsewhere. Finally we want to comment on the 
computation of the thermal pressure. It is well-known that the computation requires the 
correct normalisation of the path integral, only physical degrees of freedom contribute to the 
pressure. In the present formalism this is the correct choice of Ck,T- This amounts to the 
projection of the field independent part of the flow in momentum space onto the transversal 
parts of momenta orthogonal to the gauge fixing vector n. Such a procedure immediately 
results in 

P = iN'-l)^. (53) 

The computation of higher contributions, in particular the resummation, requires the full 
one loop running of vertices and propagators. We also emphasise that (18) already provides 
the closed expression for the resummation. We hope to report on this matter in near future. 



VIII. DISCUSSION AND CONCLUSIONS 



In this paper, we have developed a thermal renormalisation group in terms of the func- 
tional Tt[4>] = ^tI'^] ~ ro[0]- The approach implements a thermal projection at every 
integration step of the Wilsonian flow. From Tt[4'] all thermal observables - including the 
thermal pressure and the thermal self-energy - can be obtained. Our flow for Ft is free of 
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infrared and ultraviolet divergences and therefore well-suited for numerical and analytical 
studies. The running zero-temperature effective action serves as a boundary condition for 
the flow. The flow for Tt[4>] can be solved with standard techniques including expansions in 
derivatives or vertex functions, and is not confined to the weakly coupled domain. Moreover, 
the thermal flow derives from a path integral representation of the theory. This ensures that 
no double counting occurs in the flow for fT[0] and systematic solutions thereof. While 
the present construction is based on the imaginary time formalism, it is straightforward to 
implement these ideas even in a real-time formulation. 

For scalar theories, we have shown that the leading-order resummed perturbation theory 
follows independently of the scheme, and independently of the truncation. Quantitatively, 
this has been worked out to leading order within a thermal derivative expansion. Further- 
more, an optimised thermal flow for the effective potential has been provided. Its simple 
analytic form allows for stable numerical integrations, in particular at strong coupling. As 
an aside, we point out how [38] - where an infrared regulated resummation for the thermal 
pressure of scalar theories has been proposed - is related to our work. If we restrict the full 
thermal flow (12) to scalar fields, and to leading order in the derivative expansion of the 
effective action at vanishing field, then (12) corresponds to the imaginary-time analogue of 
[38]. Hence, [38] has the interpretation of a leading-order Wilsonian fiow with mass term 
cutoff. Since (12) defines the field-dependent effective action, it allows for an extension of 
[38] to higher order operators as well as to fermions and gauge fields. 

For gauge theories, the following picture has emerged: the standard Wilsonian fiow 
equation (6), equipped with the modified Ward identity (49), allows for the consistent 
computation of infrared quantities starting with an initial action in the ultraviolet. Gauge 
invariance is ensured in the physical limit. A gauge invariant implementation for all scales 
A; - in the Wilsonian sense - is problematic due to the poor performance of a mass-term 
regulator for quantum fiuctuations in the ultraviolet. The important new result is that the 
thermal fiow (10), instead, stays well-defined even for a mass-like regulator (12) since the 
zero temperature contribution renders the fiow finite. A mass term regulator fully pays off 
when combined with the axial gauge fixing, as it implies a gauge-invariant thermal fiow 
for all scales k. The difference to the flow (6) stems now from the initial condition, which 
is no longer the bare action, but the running T — quantum effective action, or some 
approximation to it. 
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